Smoking guns of a bounce in modified theories of gravity through the spectrum of the 

gravitational waves. 
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We present an inflationary model preceded by a bounce in a metric theory a la f{R) where R 
is the scalar curvature of the space-time. The model is asymptotically de Sitter such that the 
gravitational action tends asymptotically to a Hilbert-Einstein action, therefore modified gravity 
affects only the early stages of the universe. We then analyse the spectrum of the gravitational 
waves through the method of the Bogoliubov coefficients by two means: taking into account the 
gravitational perturbations due to the modified gravitational action in the f{R) setup and by simply 
considering those perturbations inherent to the standard Hilbert-Einstein action. We show that there 
are distinctive (oscillatory) signals on the spectrum for very low frequencies; i.e. corresponding to 
modes that are currently entering the horizon. 
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I. INTRODUCTION 

The inflationary paradigm while being the most ac- 
cepted and successful theoretical framework today to ex- 
plain the homogeneity, flatness and horizon problems 
docs not give any clues as to what mechanism is respon- 
sible for the origin of the inflationary era of the universe 
As such, many different model have been developed, 
from scalar fields to modified theories of gravity that try 
to describe the early accelerated expansion of the uni- 
verse that leads to the seeds of the universe we observe 
today (cf. for example In order to gain some in- 

sight on the dynamics of the early universe and to begin 
to discriminate between these models, one needs to probe 
the spectrum of the cosmological perturbations: for ex- 
ample the scalar or tensor perturbations, that arise from 
the quantum vacuum fiuctuations during the inflationary 
period. In fact since the evolution of these perturbations 
depends on the characteristics of the mechanism causing 
the acceleration of the universe, it is reasonable to assume 
that a given inflationary model will have some character- 
istic imprints on the energy spectrum of the Gravitational 
Waves (GW). 

Another issue is that the inflationary paradigm docs 
not solve the big bang singularity, in fact it is expected 
that such a singularity would be resolved within the 
paradigm of quantum cosmology @. A plausible can- 
didate regarding this is loop quantum cosmology, which 
predicts a bounce instead of the big bang singularity . 
Bouncing cosmologies have a long history and can appear 
for example in modified theories of gravity [8l-[l0|. 

In this paper, we propose a bouncing scenario within 
the context of modified theories in the metric formalism. 
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more precisely we will consider an f{R) model. These 
kind of models have been quite popular recently as they 
could describe the current acceleration of the universe 
without invoking a dark energy component, i.e. the mod- 
ification in the modified Friedmann equation mimics the 
effects of a dark energy component with the advantage 
that the acceleration is a pure gravitational effect and 
is not due to an unkown dark energy component (see for 
example |T^ - |20j ). We would like to point out that bounc- 
ing cosmologies within a Palatini formalism of /(i?) the- 
ories and connection with an effective approach to loop 
quantum cosmology have been analysed in Ref. . On 
the other hand, we would like to remind that one of the 
pioneer inflationary models was formulated within the 
context of modifled gravity |2l| . 

The bounce wc will consider in our model is followed 
by an inflationary era which is asymptotically de Sitter 
where, in addition, the gravitational action approaches 
the Hilbert-Einstein action on that regime, such that the 
modiflcation to Einstein General Relativity (GR) afects 
exclusively the very early universe, around the bounce 
and a few efolds after that. We will constrain the model 
obtaining the spectrum of the stochastic gravitational 
fossil as would be measured today. Such an analysis willl 
be carried out using a generalization for /(i?)-gravity of 
the method of the continuous Bogoliubov coeflicients flrst 
introduced in (22i - i24) . and later applied in [25|-i27j. One 
of these Bogoliubov coefficients gives the density of gravi- 
tons of the universe. An altervative approach to obtain 
the spectrum of the GW in modifled theories of gravity 
was developed in jH-fsij 

The outline of the paper is as follows. In section |lTl 
we deflne the behaviour of the bounce through its scale 
factor and derive the appropriate f{R) action compatible 
with such an evolution of the early universe. In section 
mil we summarise the methodology used to obtain the 
spectrum of the stochastic GW which is based on the 
Bogoliubov coefficients. We present in section IIVI the 
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spectrum of the GW for the model introduced in section 
In] using both a GR setup and an f{R) setup to analyse 
the evolution of the gravitational perturbations. Finally, 
in section |V] we present our conclusions. 



II. A BOUNCE IN /(R)-METRIC GRAVITY 

In this paper we introduce a phcnomenological bounce 
in the early universe as a way to avoid the big bang sin- 
gularity. A bounce at time tf, can be characterized in 
terms of the scale factor by: 



a(tb) = and a{th) > 0. 



(2.1) 



Inspired on the dc Sitter solution for a closed 
Friedmann-Lemaitre-Robertson- Walker (FLRW) space- 
time, we define the scale factor around the bounce as: 



a(t) = ab cosh [Hinf{t - 4)] 



(2.2) 



where a{t) is the scale factor, at is a constant quantifying 
the size of the universe at the bounce. In addition, Hmf 
is related to the energy scale of inflation just after the 
bounce 0. For simplicity, we will set ti, to be zero. 

While the previous solution (|2.2p corresponds to a dc 
Sitter space-time, with spherical spatial section in GR, 
here we are looking for a modified theory of gravity of 
the kind f{R) such that Eq. (|2.2|) is a solution of the 
modified Friedmann equation of a FLRW universe with 
flat spatial sections. This is a simple and straightforward 
way of modelling a bounce within the paradigm of mod- 
ified theories of gravity. In addition, our choice is not 
only based on criterion of simplicity but also inspired on 
the fact that the above mentioned solution is asymptoti- 
cally de Sitter in the past and the future (for space-times 
with flat spatial sections). Consequently, this procedure 
assures avoiding the initial singularity and guarantees a 
period of inflation leading to the first seeds of the struc- 
tures we see nowadays. It is worthy to stress that the 
existence of such a bounce is prohibited for a flat space- 
time in GR unless exotic matter, which violates the null 
energy conditions, is invoked (32| . In this paper, we work 
within the framework of /(i?) theories of gravity which 
allow for a rich description of the universe (see for ex- 
ample [l^ - |20j ) and as we will show this approach allows 
for the kind of bounces described by Eq. (|2.2[) . For alter- 
native approaches on bouncing cosmologies in modified 
theories of gravity cf. for example Refs. |8l-[l0l|. 

We remind that the action of an f{R) gravity reads 
MM: 



S 



1 

2^ 



(2.3) 



^ Notice at this regard that a > at all times and in particular for 
t > t^. In addition, d > for t > th- 



which leads to the modified Friedmann and Raychaud- 
hury equations [H-Ill [l^l : 




+ = — [p 

JR 



2;^^ -^"^ 



fjR) - JrR 
2n? 



(2.4) 



!rRR^ + !rrR + '^H fiiB.R 



(2.5) 



respectively, where H is the Hubble parameter, = 
SttG; G is the gravitational constant, p and p arc the 
energy density and pressure of the matter content of the 
universe. The dot stands for the derivative with respect 
to the cosmic time and the R subscript indicates a deriva- 
tive with respect to the scalar curvature. 

Using the above written equations (see Eq. (|2.4p ). we 
can deduce the form of f{R) such that the scale factor 
scales with the cosmic time as shown in Eq. (|2.2p . More 
precisely, from Eq. (|2.2p wc derive the temporal evolution 
for the Hubble parameter, iJ, and the scalar curvature, 
R, in a spatially flat universe: 

i/(t) =i7inftanh(jyi„ft), (2.6) 



R{t) = 6Hli [1 + tanh2(i/i„ft)] . (2.7) 

Wc will assume that the inflationary era is induced 
by modifications of gravity with respect to GR. This is 
a quite old idea first proposed by Starobinsky in its in- 
flationary model [2l|. Therefore, we will assume that 
the energy density, p, and the pressure, p, are zero. Fi- 
nally, by inserting Eqs.([2?2|. (p^ and ([^J]) in the mod- 
ified Friedmann equation (j2.4[) . wc obtain a constraint 
on f{R), more precisely / must satisfy a second order 
differential equation: 



Hf + {H^ - 2H)f + 2HHf = 0. 



(2.8) 



It is worthy to stress that this equation has two solu- 
tions [l^ - |20j . Therefore, there is not a unique form of 
f{R) given a specific evolution for the scale factor. We 
will have to impose some physical criteria to choose the 
appropriate one. The solution of equation (j2.8p can be 
written as a linear combination of the functions fi (t) and 
/2(t): 



V3, 



hit) = V6tanh^(-ffinfi) + 3cos + ^2(0 



/2(i) - y^6tanh2(i/i„ft) + 3sin (^1^61(1) + 62(1) j , 

(2.9) 
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where 9i {t) and 82 {t) are given by: 

6*1 (t) = sgn(t) arccos(l - 2ta.TLh^{H,^tt)) 



2it) 



arccos 



/ tanh^(iJinft) 
2tanh^(iJinft) + 1 



(2.10) 
1 



Here the sgn function is defined by sgn(t) = 2u{t) 
where u{t) is the Heavysidc step function |33| . 

We will choose the linear combination of /i and /2 
such that the effective gravitational coupling G^'^^^^ is 
positive, or equivalently fR>0 [G^'^^^'^ = G/ fu, as can 
be easily seen from (|2.4p ) . In order to impose this condi- 
tion it is easier to rewrite /i and /2 in terms of the scalar 



curvaturfl 



2771): 



V3, 



MR) = jR/Hl,-'Scos :-0,{R) + e2iR) 



R/Hl,- 3 sin I ^9,iR) + e2{R)] • (2.11) 



(2.12) 



/2(i) = 

Therefore, the solution for f{R) is given by: 
f{R)^C,,h{R) + C2f2{R)- 



Differentiating (|2.1ip and imposing //^ > we obtain 
the following conditions for the linear coefHcients Ci and 
C2: 




Ci < and ^ = 
Gi 

We pick Ci and G2 to be of the form: 
and G2 



(2.13) 



r r M 

Ci ~ G cos I 



Csinf^Tr), (2.14) 



with C a positive constant. It can easily be checked that 
this choice of coefficients satisfies Eq. (|2.13p and leads to 
the following expression for the f{R) solution: 



f{R) ^GJR/H. 



X cos 




arccos 



2R/HI,- 



(2.15) 



We are left with a single constant G which can be fixed 
by imposing that asymptotically, well inside the infla- 
tionary era: (i) we recover the Hilbcrt-Einstein action 



^ Notice that wc can invert Eq. (2.7) analytically. For simplicity, 
we did not include the inverted relation between the cosmic time 
and the scalar curvature. 



and (ii) the gravitational coupling reduces to the gravi- 
tational constant G. Following this procedure we obtain 
G = 2H^^f. In Fig. [U we plot the solution (|2?T5l) and its 
first two derivatives with respect to the scalar curvature. 
Well inside the inflationary era, R « 12H?^f, the sec- 
ond and higher derivatives of f{R) become small and the 
action can be approximated as f{R) « i? — 6i??jf, there- 
fore we recover the standard behaviour of GR. However, 
for low values of the scalar curvature, R, the deviations 
from GR become more significant, in particular close to 
the bounce. In fact at the bounce the second derivative 
and higher derivatives of /(i?) get very large. Despite 
the divergence of fjin, and therefore of fB.RB., the Fried- 
man and Raychandhuri equations, Eqs. ()2.4p and ()2.5p . 
are well defined for all R. We will next show, however, 
that this behaviour does have serious implications espe- 
cially in the scalar sector. 

It is well known that metric f{R) gravity can be re- 
cast as a scalar-tensor theory, in particular a Brans-Dicke 
theory with Brans-Dick^l parameter wq = and a new 
scalar degree of freedom, the scalaron, defined as (f) = fa 

14 , l33-[39ll . The squared mass of this field is given by 

38,Efl: 



Jr. - R/rr 



3/ 



RR 



(2.16) 



The squared mass of the scalaron in our case vanishes 
asymptotically, roughly when the scalar curvature is 
about R « 8.72i/inf. Unfortunately, it becomes negative 
for smaller values 



III. ENERGY SPECTRUM OF 
GRAVITATIONAL WAVES 

We now analyse the gravitational wave (GW) spec- 
trum predicted by the model introduced in the previ- 
ous section. In order to obtain the current energy spec- 
trum of gravitational waves, we use the generalization 
for /(i?)-gravity of the method of the continuous Bogoli- 
ubov coefficients first introduced by Parker, Ref. [22|, and 
Starobinsky, Ref. [2^ . Here a Bogoliubov transformation 
is applied to quantify the tensor perturbations of the 
metric in terms of time-fixed annihilation and creation 
operators and their linear coefficients a and /3. These 



^ Palatini f{R) gravity is also a Brans-Dicke theory with Brans- 
Dicke parameter ojq = 3/2.ll(]|. l38ll 

* It is usually assumed that the condition that the scalaron is not 
a ghost or a tachyon is closely related to the Dolgov-Kawasaki 
instability [T^. l4lj. which implies fun > 0. Notice that in our 
case despite having fjiji > 0, the squared mass of the scalaron 
becomes negative around the bounce. The reason behind this 
behaviour is that close to the bounce the model deviates strongly 
from GR, while in deducing the condition /rh > to avoid the 
Dolgov-Kawasaki instability, it is assumed a small deviation from 

GR MM- 
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R/^^inf R/^inf R/Wjnf 

(a) (b) (c) 

FIG. 1. These plots show from left to right: (a) the behaviour of f{R) as a function of R/H?^f (see the blue curve) and the 
Hilbert-Einstein action R — 6H^^{ (see the red curve) corresponding to the limiting behaviour of f{R); (b) the behaviour of fn 
as a function of R/Hf^f, (c) the behaviour of fnn as a function of R/Hi^f. 



coefficients, which carry the time dependence of the vac- 
uum state, fulfill the initial conditions a(tini) = 1 and 
/3(iini) = 0. Furthermore it can be shovifn that gives 
the density of gravitons of the universe. A second linear 
transformation is introduced to describe the evolution of 
the graviton density b y th e more practical X and Y vari- 
ables, which obey (IlHlllil: 



X' = -ikY, 



Y' 



a 
a 



X. 



(3.1) 



(3.2) 



Here a prime indicates a derivative with respect the con- 
formal time, r] (a = dt/drj). These new variables are 
related to the graviton density by |/3p = |X — and 
must respect the conditions |X(tini)| = |i^(iini)| = 1 and 
Re{X- Y) = 1. We highlight that in a flat FLRW de Sit- 
ter space-time, the system of equations p.ip and p.2|) 
admits the solution 1271: 



Using the definition of z we can expand into a" /a 

plus a correction term due to the f{R) modification, S, 
that vanishes at the end of the f{R) era: 



z a a fn 2 fii 
a 



Ik 

fR 



(3.7) 



Note that even though the definition of S contains terms 
that diverge at the bounce because they involve fuR and 
higher derivative^, S is always finite. 

The definition of the scale factor during the early uni- 
verse, See Eq. (|2.2I) . as well as the solution obtained for 
f{R), Eq. ([233, allows us to cast both terms of Eq. ((3J)) 
as functions of the cosmological time: 



-{t)=Hlfal [2cosh2(iJi„fi)-l] , 



(3.8) 



Xia)= l + ^ — 



.aH\ ,_k_ 



Yia) 



(3.3) 



(3.4) 



The graviton density, at the present time, 770 , de- 
fines the dimensionlcss relative logarithmic energy spec- 
trum of GW for the frequency uj as [l^. El] : 



^GwiuJ,rio) 



^'1/3(^0)1 = 



(3.5) 



where c is the speed of light in vacuum and h is the 
reduced Planck constant. 

In /(i?)-gravity the same method can be applied by 
replacii ig a" /a in Eq. p.2p by z" / z where z ~ (1\/Jr 

[IlllliHiili: 



Y' ^-^ik^ --]X. 



(3.6) 



E{t) =i/i>^ j4 - 5cosh"(i?i„fi) 
-I- [9coah^{Hinit) - 6] Z{t) 



[3cosh\Hin{t) - 2]' 
4sinh2(iJintO 

where the Z{t) is defined as: 



■Z\t) 



3 sinh (Hinit) 



cot 



(3.9) 



(^01 (t) + 02(t)) + ^J3sm\l^{Hin^t) 



(3.10) 



The function Z(t) is a monotonic increasing function of \t\ 
that goes from Z{0) =0 to the limiting value during the 



^ It can be easily checked that and are proportional to fuR 
and /rrr. 
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FIG. 2. This Fig. shows: (i) the potential z" j z (see the 
blue curve) and its asymptotic behaviour (see the blue dashed 
line); (ii) the potential a' ja (see the purple curve); (iii) the 
comoving wave- number a^H^{see the red curve); as functions 
of cosmic time and near the bounce. The variables are pre- 
sented in Planck units 



FIG. 3. This Figure shows: (i) the potential z" /z (see the 
the blue curve); (ii) the comoving wave-number a^H^ (see 
the red curve); as functions of the scale factor. The variables 
are presented in Planck units. 



dc Sitter like expansion Z{t — > oo) = 2/3. In this regime 
S « H?^fa^/S, therefore, we can make the approximation: 



■(^»^i;f) 



Hlfal [2 cosh2(iJinft) - 2/3] . (3.11) 



universe throughout the different eras as follows: 



[A 



l-(^) 

B 



4(1 + q) J 

4 

+ Pn 



l + Q 



4pA 



f{R) era, 

mGCG era, (3-13) 
ACDM era. 



In Fig. [21 we plot the behaviour of a" /a and z" /z near 
the bounce. The striking feature of the f{R) modifica- 
tion is the fact that z" /z becomes negative very close to 
the bounce {z" /z < approximately for \t\ < 0.22Hr^^). 
We note that, in GR, the potential a" /a can become neg- 
ative whenever p — 3p < 0, for example for stiff matter. 
The deviations from the GR term are maximum at the 
bounce: max|S| = |S(0)| = A.7A278H^^fal. 

When the universe is well inside the inflationary era 
the deviations of the action of our model from a Hilbert- 
Einstein action become small and equations (j2.4[) and 
(|2.5p approach the form of the Friedmann and Raychaud- 
hury equations derived from GR with an effective cos- 
mological constant. We make use of this asymptotic be- 
haviour by switching entirely to a GR description at this 
point which allows us to describe the late time evolution 
of the universe by the well known ACDM model. To ob- 
tain a smooth transition between the ACDM model and 
the model introduced in Sec. [ill we use a modified Gener- 
alized Chaplygin Gas (mGCG), see for example Ref. (isj 

(also Hi-ia). 

The effective energy density in vacuum during the f{R) 
era can be expressed as: 



2K^fR 

3iJ2 



RR 

Jr 



R 



(3.12) 



Piecing it together with the densities for the mGCG and 
ACDM models we can express the energy density of the 



The constants A and B are fixed by comparing the 
asymptotic behaviour of the density in mGCG with these 
in the f{R) and ACDM regimes: 



. l+a 



B = {pra^} 



(3.14) 



In the later mGCG and ACDM eras, the correction term, 
S, vanishes and z" / z is reduced to the GR term which is 
determined as: 



— = -^a (P- 3p). 
a 6 



(3.15) 



Here the energy density is defined by the second and 
third lines in the definition p.l3p . while the pressure is 
determined by the equations of state of the respective 
models: p = (p - iA/p")/^ for the mGCG; p = T,ip, = 
T.WiPi for the ACDM model where w = 1/3 for radiation, 
w = for dust and w = —1 for the cosmological constant. 

Combining these results to write dowro z" / z during 
the three different eras we obtain: 



Hi, 



2a2 -al) +E 



z 
z 



K^o?A {A 



, a < ai, 

, a\ < a < ail 

, a > 02. 



(3.16) 



The transition from the first to the second era takes place 
at a = ai, and from the second to the third era at a = 02. 



Notice that z" j z = a!' ja for a > ai 
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The values of both ai and 02 are imphcitly defined by the 
condition that the branches of Eqs p.l6p are continuous 
at the transition time. 

To determine the gravitational energy spectrum flaw 
we integrate equations (j3.ip and (j3.6p for a given wave- 
number k from an initial time tmi set before the bounce 
until the present time. The value of the parameter 
iini defines how much the universe contracts before the 
bounce. The wave-numbers k are scanned from the mode 
that is entering the horizon at present, fchor = 0,0 Hq, 
to the mode that corresponds to the maximum of z" /z, 
kmax = ■\/max(z"/z). In order to understand the im- 
portance of the f{R) corrections in equation p.6p we 
determine the evolution of X and Y using two methods. 
In the first method the integration is made considering 
solely the GR term a" /a while in the second method the 
complete definition of z" /z is used. 

In order to solve the system of differential equations, 
Eqs. p.ip and (|3.6|f 1. we need to set proper initial condi- 
tions. Even though the complex expression of the f{R) 
term z" / z gives little hope of obtaining an exact ana- 
lytical solution, we can show that, when only the GR 
term, a" /a, or the approximation during the de Sitter 
like regime are considered, the equations can be solved 
analytically. This allows us to obtain initial conditions 
for the numerical integration which we will apply on the 
next section. 

We begin by combining the two first order differential 
equations, p.ip and p.6p . that govern the evolution of 
the variables X and y in a single second order differential 
equation for X: 



X' 



e- — \x. 

z 



(3.17) 



Notice that Eqs. p.Sp and p. lip can be cast in the gen- 
eral form gy. 



(3.18) 



with 7 a parameter taking the value 7 = 1 for the GR 
case (gi = a" /a, c.f. Eq. ([3^) and 7 = 2/3 for the 
limiting de Sitter-like regime within the f(R) framework 
(52/3 = z"/z{t > Hrl), c.f. Eq. (|XTT]) ). Substituting 
g-y in p.l7p and expressing the differential equation in 
terms of the cosmic time we obtain: 

X+Hinfta-nh{Hinft)X 

+ Hi, [2 cosh2(i7i„ft) _ 7 _ ^ = 0, (3.19) 

where the reduced wave-number q is defined as g = 
k/{abHini). Equation p.l9p has a solution given by the 
linear combination: 



X{t) = DiXi{t)+D2X2{t), 



(3.20) 



^ Notice that z" / z reduces to a" /a when H vanishes. Therefore, 
we can unify in Eqs. 1 13.11 1 and 113.61 1 the approaches we will follow 
in the GR and the f{R) frameworks. 



of the real- valued functions Xi and X2: 



X cos 



X sm 



^sinh^ {Hinit) 


+ 7 + 92 


7 + 172 




gn(7 1 + g^) 


i(a)\/7- 


^sinh' (Hinit) 


+ 7 + 92 



7 + r 



sgn(7 -l + q'^) 



(3.21) 



The phases 01 and 02 are defined as: 

4>i{t) = TT + sgn(t) arccos [2sech^(i/inf<) — l] , 

2(7 + 9') 



) (i) = TT + sgn(t) arccos 



sinh^(i7i„ft)+7 + g2 



(3.22) 



The corresponding Y functions can be obtained by in- 
serting Eq. (f33T|) in Eq. (|3?T|) : 



Yiit) = - 



- V7 + 



sinh(i7inft) cosh' (i/inft) 
sinh2(iJi„ft) + 7 + ^2 

|7-l + g2| 



sinh'(ifinft) + 7 + 9' 
sinh(_ffi,ift) cosh' (i/inft) 



+ V7 + 



sinh'(iJinft) + 7 + 9' 
|7-l + g2| 



sinh^(FinfO + 7 + 9' 



Xi{t) 



(3.23) 



The linear coefficients Di and D2 are determined by 
the initial conditions |X(iini)| = |F(iini)| = 1 and by im- 
posing that, for high fc, X converges to a Bunch-Davies- 
like solution [3| with negative Hubble parameter, during 
the contraction phase. 



(3.24) 



Do = 



V -''^l (^ini) + ^2'(^ini) 

q XsiUnOXiiu^o + ri(ti„i)i2(tini) 



\l - I + q^lV^iT^ 



\/ ^'(iini) + Y2(t\n\) 



(3.25) 



By fixing the initial conditions before the bounce, we 
are effectively defining a range of wave-numbers, (fc2 < 
2:"/z(tini), that at the beginning of the integration are 
already outside the horizon. For those modes there is an 
increase in the graviton density, |/3|2, before the bounce. 
This fact is reflected in an upward shift in the energy 
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spectrum, flaw, ^-s we will show. As z"/z grows expo- 
nentially with negative time, the further away from the 
bounce that we set the initial conditions, the more intense 
this shift is expected to be and the more wave-numbers 
are affected. 



IV. NUMERICAL SIMULATIONS 

In our model we can identify a set of four independent 
parameters: 

1. the time when the initial conditions are set, ti^i 
which fixes the relation Oini/ab between the initial 
value of the scale factor, Oini, and the scale factor 
at the time of the bounce, ab; 

2. the scale factor at the time of the bounce, Ub; 

3. the energy scale during the inflationary era, i?inf, 
which is related to the constant iJinf by: 



inf ■ 



k2 '"f 



1/4 



(4.1) 



4. the value of the parameter a in the mGCG model. 
As this parameter changes the spectrum only in 
the high frequency range, as opposed to the effects 
of the bounce which are relevant only for the low 
frequencies, we use the value a = —1.04 obtained 
in Ref. [i^ which corresponds to the best fit using 
the WMAP data. 

For a chosen set of values {ab,tini, Einf) we compute 
the GW spectrum via the numerical integration of the 
differential equations p.ip and p.6p in two ways: case 
(i) within a GR perturbative treatment; case (ii) within 
a full f{R) treatment. The integration is divided in two 
parts with the first one, during the f{R) era, performed 
in terms of the cosmic time and the second one, during 
the mGCG and ACDM eras, performed in terms of the 
scale factor. Following previous works, [13,IH,S,[131j 
stop the integration when 3> z" /z, i.e. when the mode 
is well inside the horizon, since in this regime X and Y 
have a highly oscillatory sinusoidal behaviour while 
is virtually constant. This is an excellent way of reducing 
the computing time of the spectrum. 

Before describing the two methods we define two par- 
ticular wave-numbers. The first one: 



kj — ^ymax[^\ = 



(4.2) 



can be interpreted as the smallest wave number for which 
the solutions in Eqs. ((X^ and ([5:^ (with 7 = 2/3) are 
a good approximation during the entire f{R) era. The 
second one: 



kii = \/z"/zitini), 



(4.3) 



indicates the mode with higher wave number that sees its 
graviton density increased during the contraction phase. 



Notice that, in principle, all modes should have an in- 
crease on its graviton number before the bounce and, 
indeed, if one sets tmi further away from the moment of 
the bounce the limit fc// grows. However, by starting 
the integration from a given initial fixed time for all the 
modes only those modes close to the bounce will have a 
significant increases on its graviton density. 

We proceed to describe the two methods employed to 
obtain the energy spectrum of the gravitational waves: 

• In case (i) only the GR term a" /a is considered 
when performing the integration, i.e. we integrate 
the E qs. ([XT]) and ([3^. Since Eqs. ((33T|) and 
(|3.23p (7 = 1) provide us with analytical solutions 
during the f{R) era, we can skip the numerical in- 
tegration all together during this period and use 
Eqs. p.2ip and p.23p to set the initial conditions 
for the integration during the mGCG era. This is 
done in two different ways depending on whether 
the mode has already crossed the horizon at the 
moment of the transition defined by ai or not (we 
remind that ai defines the transition from the f{R) 
solution to the mGCG era (cf. Ea. (|XTC)) ). 

— If the mode k is outside of the horizon at 
the beginning of the mGCG era {k < aiHi), 
we use Eqs. (|3.2ip and (|3.23p to fix the ini- 
tial conditions for the integration at a; = ai 
and then integrate numerically until the mode 
reenters the horizon, more precisely until a/ ~ 



lO^Cg^tj-y, where al'^^j.y is the value of the scale 
factor at the moment of the reentry of the 



mode defined as k = a';,'^^^^!! (al'^^^.y) . As said 
previously, we stop our intregration slightly af- 
ter the mode reenters the horizon because, for 

(k) 

a > a^ntryy ^iid Y have an oscillatory be- 
haviour; thus, \f3\ remains constant after the 
crossing. 

For higher modes {k > aiHi), which are 
still inside the Hubble radius at the moment 
of the transition, we consider that the solu- 
tions (|3.2ip and (|3.23p are valid slightly be- 
fore the mode crosses the horizon, i.e., until 

Gi 10~^a^^|^ , where a|,^jj. is the value of the 
scale factor at the moment when the mode ex- 
its the horizon defined by fc = aitit^('^cxit)- 
is noteworthy to point out that, for large wave 
numbers, the solutions obtained for this model 
are virtually identical to the dc Sitter solutions 
after the bounce and this is precisely the rea- 
son why we can apply them during the mGCG 
era when dominated by the effective cosmolog- 
ical constant A. Therefore, we use these solu- 
tions to set the initial conditions at a = and 
perform the integration until a/ ~ 10^a^Jjjj.y, 
when the mode is again well inside the hori- 
zon. The reason we do not impose for the 



(fc) 



8 



modes aiHi < k) the same boundary condi- 
tion we used previously for k < aiHi) is to 
avoid the oscilatory regime before the modes 
cross the horizon as it is very time consuming. 

• In case (ii) the fuh definition of z" / z is considered 
(see Eqs. (|3.7p and p.9p ). therefore we integrate 
Eqs. (|3.ip and p.6p . The integration is performed 
accordingly to one of the following procedures: 

— If (fc < fliiJi) and [k < kj), we begin the in- 
tegration during the f{R) era by setting the 
initial conditions at t = ti^i using Eqs. p.2ip 
and ([3:231 (7 = 2/3). The system of differ- 
ential equations is then integrated until the 
moment of the transition at t — ti, where ti 
is defined by ai = at cosh(i7infii), as given in 
Eq. (|2.2p . We then use these results as ini- 
tial conditions for the integration during the 
mGCG era and then integrate until after the 
mode reenters the horizon at ^ 10^a^^|.jy. 

— If (fc < ayHi) and (fc > fc/) the approxi- 
mate solutions given by Eqs. p.2ip and p.23p 
(7 = 2/3) are not affected by the exact form of 
the potential z" / z during the bounce. There- 
fore, we can skip the numerical integration 
during the /(i?) era, using those solutions to 
set the initial conditions for the integration 
during the mGCG era at ai = oi and then we 
integrate numerically until the mode is well 
inside the horizon at a/ ^ 10^a^^|.j.y. 

^ If (^ ^ o-iHi) and (fc < fc/) we again set 
the initial conditions for the integration of 
the f{R) period at t = tini using Eqs. p.2ip 
and (|3.23p (7 = 2/3). The numerical in- 
tegration is then performed in terms of the 
cosmological time until the mode approaches 
its horizon exit, i.e., t ~ ti defined by ai = 
Of, cosh(iJinfii) ~ 10~^a|,^jj. These results are 
then used as initial conditions at a = and 
the integration is continued for the mGCG 
and ACDM eras the mode are well inside the 
horizon at a/ ~ 10^0^^^^.^. 

^ If (^ Jc^ o-iHi) and (fc > fc/) we can use the 
approximated solutions Eqs. p.2ip and p.23p 
(with 7 = 2/3) to set the initial conditions 
for the integration of the mGCG era before 

the mode exits the horizon at ~ 10~^agxjf 
The integration is then performed until a/ ~ 



2„(fe) 
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when the mode is again well inside 



the horizon. 

The range of frequencies scanned is delimited by 
Wmm = Koi/aQ « 1.43 X 10"^'' rad/s and LOmax = 
kmax/cLo « 1.49 X 10~-^^£'inf rad/s. The values of the 
parameters of the ACDM model were taken from the 
WMAP7 data [U: fir = 8 x 10-^ = 0.272, 



= 0.728 and Hq =70.4 km/s/Mpc. The value of 
scale factor at present time is set to oq = 10^®. 

Finally the following constraints on the energy spec- 
trum of GW were considered [i^, [s^] : 

• From the CMB radiation: h^flGwi^horiVo) ^ 7 x 
IQ-" for ho = i?o/(100 km/s/Mpc) and ujhor = 
2 X lO-i^/io rad/s; 

• From observation of milliseconds pulsar: 
hlflcwi^^puhm) < 2xl0~^ forwpij =2.5xl0"^/io 
rad/s; 

• From the Cassini spacecraft: hQflGwi^Cas^Vo) < 
0.014 for ujcas = 7.5 x 10~^ho rad/s; 

• from the LIGO experiment: hQilcwi^, Vo) < 3.4 x 
10^^ for frequencies on the order of a few hundred 
rad7s; 

• From BBN: hln{Lj,To)duj/uj < 5.6 x 10"^ for a;„ w 
lO-^rad/s. 

The results obtained for the GW spectra are shown 
in Figs, m [5] and |6l In each figure two of the three pa- 
rameters (aini/ab, Ob, -Einf) are fixed while the third one 
is changed so as to illustrate its effects on the spectrum. 
Furthermore, in Fig.|4]the results obtained using the per- 
turbative approach in (i) GR and (ii) f{R) are compared. 

As expected, the existence of a bounce in the early 
universe affects the spectrum only in the low frequency 
range, where a highly oscillatory regime is present in 
contrast with the smooth plateau in the intermediate 
frequencies and the rapid decay in the high frequency 
range. The fact that the oscillatory structure appears in 
the spectra of both treatments (i) and (ii) suggests that 
it is due to the existence of the bounce and not a conse- 
quence of the effects of /(i?)-gravity. Similar oscillations 
have been obtained in works of loop quantum cosmology 
first pointed out by Afonso et al [531, |5J] . 

The differences in the spectra obtained in (i) a GR 
setup and in (ii) an f{R) approach are evidentiated in 
Fig. I4b[ These appear only for very low values of the 
reduced wave-number, q, which corresponds to the low 
frequencies of the spectra where an extra local minimum 
occurs for fc = ObiJinf/v^ . This corresponds to the value 
that makes the factor sgn(7 — 1-1- q^) = in solutions 
([5:^ and ([X^ for 7 = 2/3. Notice that for the medium 
and high frequencies, and even for the higher frequencies 
on the oscillatory regime, the spectra obtained with the 
two treatments are virtually indistinguishable. The spec- 
tra obtained are in conformity with the assumption that 
the maximum value of the f{R) corrections, S, sets the 
scale on the upper limit of the wave-numbers for which 
there are noticeable differences in the spectra determined 
using the methods (i) or (ii). Thus, for the gravitational 
spectrum flaw to be useful in discriminating between 
the GR and f{R) theories describing the evolution of the 
tensor perturbations, the wave-number fc/ must be within 
the accessible range of wave-numbers today., i.e.: 



fc/ > fcho 



(4.4) 
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FIG. 4. Spectra of the energy density of GW at present time for different values of the scale factor at the bounce, at- The 
spectra are plotted for (a) the entire range of frequencies and (b) the low frequency range. A comparison is made between the 
results obtained using the GR setup (dashed curves) and the f{R) treatment (continuous curves). The value of ai, increases from 
the red curve to the blue curve: Ub = 2 x 10^; at = 2 x 10^; at = 2 x 10'*; at = 10 x 10^; at = W x W^. {Eini = 1.5 x 10^'^GeV. 
flini = Wab) 




w (rad/s) oj (rad/s) 

(a) (b) 



FIG. 5. Spectra of the energy density of GW at present time for different values of the scale factor at the bounce, ab- The 
spectra are plotted for (a) the entire range of frequencies and (b) the low frequency range. The value of a\n\/ab increases from 
the red curve to the blue curve: aim = lOa;,; aini = 2Qab\ Omi = SOa;,; Oini = lOOa^. (-Einf = 1-5 x lO^^GeV. at = 2000). 




FIG. 6. Spectra of the energy density of GW at present time for different values of the scale factor at the bounce, ab- The 
spectra are plotted for (a) the entire range of frequencies and (b) the low frequency range. A comparison is made between 
the results obtained with a fixed value of at (continuous curves) and a fixed value of abH^ni (continuous curves). The value 
of Ei^{ increases from the red curve to the blue curve: Ei^f = 1.5 x lO^GeV; Ei^t = 0.5 x lO^^'GeV; -Einf = 1.5 x lO^^GeV; 
Eint = 0.5 X lO^'^GeV; Ei^! = 1.5 x 10^'^GeV. (ain; = Wab. at = 2 x 10^). 

In Fig. |4] a horizontal displacement is observed in the the fact that the solutions p.2ip and (|3.23p do not de- 
spcctrum when the value of is modified. This reflects pend directly on the wave-number k but instead on the 
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10^ 10^ 10^ 10^ 10^ 
at 

FIG. 7. This picture shows the maximum allowed value of 
cbini/cLb that does not violate the CMB limit, as a function of 
ah and for different values of the energy scale, i?inf. From the 
red curve to the blue curve: i?inf = 1.5 x 10^*GeV; -Einf ~ 
0.5 X IQi'^GeV; E^^f = 1.5 x lO^^^GeV; Si„f = 0.5 x lO^'^GeV. 
(A: = 1.074a6JYinf) 



reduced wave- number q = k/{abHin{). As such it is ex- 
pected that the position of the peaks on the spectrum to 
be affected the by changes on the product ObiJinf. This 
effect is also seen in Fig. [51 because lowering the energy 
scale for a fixed value of abHinf only shifts the spectrum 
downwards without changing the position of the peaks. 

In Fig. [5] an increase on the energy spectrum of GW, 
^lGW, as well as in the range of the affected frequencies, 
is observed as the initial conditions are set further away 
from the bounce (higher Oini/af,). This results agree with 
the analysis done previously in Sec. IIIII 

The GMB constraints set an upper limit in the low 
frequency end of the spectrum of flcwii-^hor, Vo) ^ 1-4 x 
10^^". This value limits the amount of contraction that 
one takes into account as the intensity of the peaks in the 
energy spectrum scales rapidly with increasing aini/a&. 
In fact from Fig. [S] it is possible to observe that if the 
condition (|4.4p is met the relation aini/flb is heavily con- 
strained: for aini/flfc = 20 the maximum of the poten- 
tial is slightly above 10~^°, while for flini/ofc — 50 the 
maximum of the potential is almost two orders of mag- 
nitude higher which clearly violates the CMB limit. In 
order to obtain a better picture of the restrictions set 
on Oini/afc, we select the wave-number of the most left- 
ward peak (of the order of kp « 1.074af,i?inf) and re- 
peat the integration with varying aini/a;,. The results 
obtained are presented in Fig. [7] for varying energy scales 
during inflation. Although the maximum allowed value 
of Oini/ab increases slightly with a?,, it rapidly becomes 
stagnant. Notice that by lowering the energy scale, the 
limiting value of the relation aini/of, is shifted upwards, 
and with E'inf = 0.5 x lO^^GeV we are able to obtain 
Oini/af) ~ 2100. This limitation in the amount of con- 
traction considered may be avoided by lowering the prod- 
uct flbTJinf, thus shifting the higher peaks to the invisible 
region of frequencies lower than Wmin- This, however. 



means that one would lose the ability to differentiate 
between the GR and the purely f{R) imprints on the 
spectrum. 

In Fig.|n]the effects of changing the energy density dur- 
ing inflation may be observed. Besides the usual vertical 
shift on the intermediate plateau and variation of the 
maximum wave-number, the value pinf affects the posi- 
tion of the peaks on the spectrum through the constant 
o-bHirii as discussed above. 



V. CONCLUSIONS 

In this work we investigate the effects that the pres- 
ence of a bounce in the early universe has on the energy 
spectrum of the gravitational waves. In order to avoid vi- 
olations of the null energy condition that usually appear 
in bouncing FLRW cosmologies in GR, [32| , we work 
within f{R) gravity, fixing the desired behaviour for the 
scale factor and then solving the Friedmann equations to 
obtain the function / (designer f{R) gravitv |12| - [20| ). A 
mGCG model, |43l - l49j . was used to connect the bounc- 
ing and subsequent inflationary era with the the radiation 
epoch as it provides an easy way to describe a smooth 
transition between the two phases (inflation and radia- 
tion) without affecting the low frequencies of the spec- 
trum. 

The method of the Bogoliubov coefficients, [12, [25l - 
\2T\ . l42j , was used to determine the evolution of the gravi- 
ton density while treating the tensorial perturbations 
both in a GR approach and in a full /(i?) treatment, 
[3 [s^ - ISTj . The energy spectrum of the gravitational 
waves was determined for different values of the scale 
factors at the time of the bounce, different amounts of 
contraction before the bounce and for different energy 
scales during the inflationary era. 

We show that the presence of the bounce creates vari- 
ous peaks in the low frequency range of the spectrum, 
whose position and intensity depend on the relations 
flfc-ffinf and aini/ofc. We remind that Ob defines the value of 
the scale factor at the time of the bounce, Hmf is related 
to the energy scale during inflation by H?^f = ^Ef^^ and 
Oini defines the value of the scale factor at a fixed time 
before the bounce when the density of gravitons is neg- 
ligable. An increase of abH-mf shifts the peaks rightward 
on the spectrum while an increase in aini/ab results in 
an enhancement of the peaks that correspond to wave- 
numbers below fc// (cf. Eq. (|4.3|) '). The intermediate and 
high frequency regions are left unaffected by the bounce. 
The condition (|4.4[) . which we have assumed, ensures that 
the modes entering currently the horizon have some im- 
prints of the f{R) era. 

The results obtained are within all the observational 
constraints considered except for the CMB radiation data 
as, for the higher energy scales even a contraction of two 
orders of magnitude before the bounce makes the spec- 
trum violate the limits imposed on the low frequencies of 
the spectrum, (42l . [s^ . An estimation on the maximum 
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amount of the allowed contraction without violating the 
CMB constraints was determined by studying the de- 
pendence of the height of the most intense peak of the 
spectrum on the parameters of the model. The fact that 
the contraction before the bounce his highly constrained 
may be viewed as an indication of the existence of an ex- 
tra state of the universe that precedes the dc Sitter like 
contraction. The potential imprint of such a hypothetical 
state on the gravitational waves spectrum is not studied 
in this work. 

Furthermore, wc note that the results obtained in this 
work are hard to be detected in the near future. Most 
of the planned experiments that might be sensible to the 
cosmological gravitational waves operate on frequencies 
of > IQ-^Hz (see Fig. 2 of Ref. [13 and Fig. 6 of 
Ref. (ssj). while the imprints originated from the bounce 
are observable, at best, up to frequencies of < lQ~^Hz. 
The best hope for the detection of these effects relies in 
future measurement of the B-mode polarization of the 
CMB radiation 



Wc would like also to point out that the solution ob- 
tained for f{R) is differentiable for all the values of R but 
its second and higher derivatives blow up at the bounce. 
Although this divergence does not have serious implica- 
tions in the tensorial perturbations, it does imply that 
the scalaron, [iJ: Is^ - IstI . l40j , becomes a tachyon near the 
bounce, so some problems arc expected when studying 
the scalar sector. We will improve this downside of the 
model in future works. 
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